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ABSTRACT 


The three elastic constants of tantalum, tungsten and molybdenum 
have been measured in the temperature range between 4.2 and 300°K. The 
1 
7 (C14 7-C) 5) for molybdenum and the 
computed bulk moduli of all three ‘elements show an anomalous temperature 


directly measured quantities C. and 


dependence. These data along with published results for other metals 
generally show a difference in temperature dependence of the bulk modulus 
between 3d and 4d-5d metals. Overlap contributions to the elastic con- 
stants are found to favor an antiferromagnetic ordering. The Debye 
characteristic temperatures calculated from elastic constants are found 


to be in substantial agreement with results obtained from calorimetry. 
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INTRODUCTION 


The elastic constants of a material determine the linear coefficients 
relating stress and strain in that material. In themselves they constitute 
an important equilibrium property. Because the equilibrium properties of 
a system are related through considerations of the system's total internal 
energy, knowledge of the elastic constants is important to a complete un- 
derstanding of thermal, electrical and magnetic behavior. 

The elastic constants of a solid are a measure of the curvature of 
the cohesive energy and their measurement yields information about the 
interatomic forces. In particular, low temperature observations are of 
importance for comparison with the results of theoretical calculations. 
Also, the low temperature elastic constants may be used to predict the 
lattice specific heat, which may then be compared with calorimetric results. 

The transition metals have always been of great interest because of 
their incomplete d shells. The Group VB and Group VIB metals, vanadium, 
chromium, niobium, molybdenum, tantalum, and tungsten, have the particular- 
ly interesting features of exhibiting a body-centered cubic structure and 
very high cohesive energies. Until recently very few measurements have 
been made on these metals, principally because of metallurgical difficulties 


in obtaining pure single crystals. 





EXPERIMENTAL PROCEDURE 


The adiabatic elastic constants of tantalum, tungsten, and molybdenum 
have been obtained from room temperature down to liquid helium temperatures 
by measuring the transit times of elastic waves in single crystals. Un- 
rectified waveforms were observed using an ultrasonic apparatus similar 
to that described pEemieiedy 2A preliminary report on the elastic con- 
Stants of tantalum pointing out anomalies in the observed bulk modulus 
has already been given 

The tantalum specimen was purchased from Materials Research Corpo- 
ration as a quarter inch diameter rod from which a 2.0643 cm (at 20°C) 
specimen was cut with accurately parallel faces. Because of its small 
diameter, no attempt was made to orient this crystal. Instead the orien- 
tation was taken "as is". This orientation was determined by Laue back- 
reflection x-ray methods to deviate from [110 | by about 8°, the asso- 
ciated polar angles being 9 = 84.9° and G = Boe7 Transit time measure- 
ments were made on the tantalum specimen both before and after annealing 
at 1200°C for five hours in a vacuum of at least ig”? mm Hg. Annealing 
did not affect the observed transit times. No change in velocity or at- 
tenuation of the elastic waves was observed in passing through the super- 
conducting transition temperature for tantalum. The transducer binder 
used with the tantalum sample was glycerin below 0°c and salol at and 
above 0°C. 

Temperatures on all samples were controlled to within 0.5°C during 
transit time measurements using a Brown Electronik recorder. A copper- 
constantan thermocouple was used for all temperature measurements above 
liquid air temperatures. For temperatures from liquid air to liquid 
helium temperatures a copper resistance thermometer was used when making 
the measurements on tantalum; a copper-gold 2.1 a/o cobalt alloy ther- 
mocouple was used for the tungsten and molybdenum measurements in the 
same low temperature range. 

The tungsten single crystal was purchased from the Linde Company 
as a 1/2" diameter rod that had been grown with its longitudinal axis 
2° off [110 | . From this a specimen was prepared that had two parallel 
faces accurately perpendicular to [110] and 2.4866 cm apart (at 0°c). 
Nonaq stopcock grease was used as the binder for tungsten at all tem- 


peratures. 
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The molybdenum single crystal was obtained from the Wah Chang Corpo- 
ration as a rough surfaced bar averaging about one inch in diameter. Using 
standard Laue back-reflection techniques a [110] direction was located 
in this specimen. Two parallel faces were then cut perpendicular to this 
direction, the mean distance between faces being 1.8322 cm at 0°c. The 
maximum deviation from parallelism of these surfaces in the molybdenum 
sample, as was also true for the tantalum and tungsten specimens, was 
0.0002 cm. The transducer binder used for molybdenum was Nonaq stopcock 
grease for the readings at and just below room temperature and Dow Corning 
High Vacuum Grease from below these temperatures down to helium temperatures. 

In calculating the elastic constants of these metals the thermal ex- 
pansion data of Nix and MeeNaiE were used down to the limits of their va- 
lidity (the liquid air temperature region). Extrapolations from these data 
following Grueneisen plots were used down to 4.2°K. The lattice constant 
used for tantalum was 3.3058A at 25 cy as given in a survey by Swanson and 
preaaee From this a density of 16.6324 grams per cubic centimeter at 0°c 
was calculated. For tungsten the lattice constant given in posreene of 
3.1650A at 25°C was used. The 0°C density of tungsten, using this value 
of lattice parameter, is 19.2710 grams per cubic centimeter. Pearson was 
also the source for the lattice constant of molybdenun, 3.14684 at 20°C. 
The 0°K density of molybdenum calculated from this value is 10.2549 grams 
per cubic centimeter. 

Results of impurity analyses conducted on these three metals are con- 
tained in Table I. The analyses of tantalum and molybdenum were conducted 
by the firm of Kennard and Drake of Los Angeles, using test material that 
had been taken from the original specimens. The tungsten analysis, pro- 
vided by the Linde Company, is an average analysis for material produced 
at the same time as the sample used. Table I shows a low metal impurity 
content for all three samples. The interstitial gas content for tungsten 


is low, larger for tantalum and quite substantial for molybdenum. 





—_——> @& @ - 
= See > <i § 6 
=e Se: «© 

_—_— ee a © _* 

- —_——_ & OS Se 
—ll aw alee eee’ & ° 
= > _— a” «—_ 

© 6 —_ =—-i'e & +<- 
—- == = i cul @i se 


=P — ee 






TABLE I. Chemical and vacuum fusion analyses of crystal specimens 
(ppm by weight) 


TANTALUM TUNGSTEN MOLYBDENUM 
Fe 20 c 10 Fe 50 
Cu 3 S 10 Si 7a 
0. o2 Al <100 0. 930 
N, 124 Be <100 N, 140 
H, 9 Cd <100 H, 10 

Mg — 

Si <100 

0. Zp 

N, < 20 

H, 1K 


All other impurities nil 





EXPERIMENTAL RESULTS 


For the tantalum specimen the measured transit times were used to 
calculate the quantity ev’ for each of the three modes being propagated 
in a direction some 8° off {110 | . Smooth lines were drawn through 
these raw data of ev versus temperature, yielding the curves of Fig. l. 
Values were taken from these curves at regular increments of temperature 
and a successive approximation method due to neeohboure: was applied to 
derive the values of the elastic constants listed in Table II. 

Since the direction of elastic wave propagation in the tungsten and 
molybdenum samples was directly along [110] » the values of ev" calcu- 
lated from the transit time measurements are equal to certain linear com- 
binations of the three adiabatic elastic constants. Accordingly, Fig. 2 
displays the measured variation with temperature for the longitudinal mode 

ev" = C. = $(C, 1 +045) + Cry, in tungsten. Fig. 3 displays the measured 
variation with temperature of the v's associated with the two trans- 
verse modes for tungsten. These equal =(C,,-C, 5) for the mode polarized 


perpendicular to the z-axis and C 4 for the mode polarized parallel to 


the z-axis. Smooth curves have ee drawn through the observed data 
points of Figs. 2 and 3 and a compilation is presented in Table III of 
the adiabatic elastic constants of tungsten for values taken at regular 
intervals of temperature from these smooth curves. 

Correspondingly Figs. 4 and 5 display the measured longitudinal and 
transverse v- for molybdenum. Table IV shows the adiabatic elastic 
constants of molybdenum derived from the curves depicted in Figs. 4 and 5. 

Taking into account the various sources of error, we estimate the 
absolute accuracy of a single observation of v- to be + 0.5 percent, 
although the internal consistency of the data is much better than this. 
Since the elastic wave propagation in the tungsten and molybdenum samples 
was directly along a principal direction, the above estimate of + 0.5 
percent should be valid as an estimate of error for the final reported 
values of the elastic constants. In the case of tantalum, however, the 
direction of propagation was some 8° off [110] . Errors attendant to 
determining the direction cosines of this direction necessitate degrading 


the estimate of absolute error of the tantalum results to + 1.0 percent. 
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Fig. 1 Measured values of ev" for tantalum as a function of temperature. 
The direction of elastic wave propagation for this sample was 8° 


removed from [110] (See Text). 
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TABLE II. The elastic constants of tantalum, in units of hg dune oaamee 
as a function of absolute temperature. Calculated by a perturbation method 
from smooth curves drawn through values obtained from measurements made in 
a direction of propagation 8° from [110] . The number of significant 
figures shows the internal consistency of the data and does not indicate 
the accuracy of the absolute value. 


es Bad op oa 5(C44-C),) Ae Cp 44 
© 2.6632 1.5816 0.8736 0.5408 1.9421 0.7080 
20 2.6633 1.5823 0.8728 0.5405 1.9426 0.7095 
40 2.6641 1.5844 0.8709 0.5398 1.9443 0.7135 
60 2.6649 1.5870 0.8682 0.5390 1.9463 0.7188 
80 2.6652 1.5896 0.8647 0.5378 1.9481 0.7249 
100 2.6638 1.5912 0.8609 0.5363 1.9487 0.7303 
120 2.6603 1.5909 0.8564 0.5347 1.9474 0.7345 
140 2.6562 1.5900 0.8522 0.5331 1.9454 0.7378 
160 2.6518 1.5888 0.8480 0.5315 1.9431 0.7408 
180 2.6467 1.5875 0.8436 0.5296 1.9406 0.7439 
200 2.6412 1.5860 0.8393 0.5276 1.9377 0.7467 
220 2.6350 1.5840 0.8352 0.5255 1.9343 0.7488 
240 2.6286 1.5818 0.8309 0.5234 1.9307 0.7509 
260 2.6221 1.5794 0.8266 0.5214 1.9270 0.7528 
280 2.6155 1.5768 0.8224 0.5194 1.9230 0.7544 
300 2.6091 1.5743 0.8182 0.5174 1.9192 0.7561 
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TABLE III. The elastic constants of tungsten, in units of lo.” dye senme: 
as a function of absolute temperature. The first three columns are readings 
taken from smooth curves drawn through the experimentally determined points. 
The number of significant figures shows the internal consistency of the data 
and does not indicate the absolute value. 


T ev). pV>. ; pY5, ; B 
i) se Cr, 3811-12? oid Cio O11 2215) ieee 
0 5.3188 1.6313 1.6380 5.3255 2.0495 3.1415 0.4182 
20 5.3188 1.6313 1.6380 §.3255 2.0495 "S.1415 0.4182 
FON) 5.3187 1.6313 1.6380 5.3254 2.0494 ammel4 0.4181 
PUM oroley 1.6312 1.6379 5.3254 2.0496 301415 0.4184 
80 5.3186 1.6312 1.6378 5.3252 2.0496 18.1415 0.4184 
[oom 5.3181 1.6310 1.6366 5.3239 2.0503 3.1415 0.4193 
120 5.3166 1.6298 1.6352 5.3220 2.0516 3.1417 0.4218 
Wo) 5.3138 1.6273 1.6315 5.3180 2.0550 3.1427 0.4277 
160 5.3100 1.6238 1.6270 5.3132 2.0592 3.1439 0.4354 
180 5.3036 1.6207 1.6220 5.3049 2.0609 1.1422 0.4402 
page 562942 1.6178 1.6172 5.2936 2.0592 3.1373 0.4414 
Page) 522646 °-1:6153 1.6124 5.2817 2.0569 3.1318 0.4413 
2200) 5.2750 1.6130 1.6077 5), 2697) 280543 aemIel 0.4406 
foome 5/2653 1.6108 1.6031 5.2576 2.0514 3.1201 0.4406 
280 5.2557. 1.6089 1.5984 5.2452 2.0484 3.1140 0.4395 
300 5.2462 1.6072 1.5937 5.2327 2.0653 jB.1078 0.4381 
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TABLE IV. The elastic constants of molybdenum, in units of ge dynes 


as a function of absolute temperature. The first three columns are readings 
taken from smooth curves drawn through the experimentally determined points. 
The number of significant figures shows the internal consistency of the data 
and does not indicate the absolute value. 





T ev, ev). evs, B 
CR) 5 Oia 3641-65) Cin Cio 7(2 4420, 9) Cron as 
gueeeaeseso 1.2503 179855 4.5002 1.7292 2.6529 0.4789 
20 4.3639 62497 1.3855 4.4997 1.7287 2.6524 0.4790 
40 4.3600 1.2488 1.3850 4.4962 1.7262 26495 0.4774 
Bam 493557 12479 1/3840 4.4918 1.7238 2.6465 0.4759 
80 4.3548 1.2468 1.3820 4.4900 1.7260 2.6473 0.4792 
iogme 443495 1.2456 iwaeo7 4.4846 1.7232 2.6437 0.4776 
fom 499389 1.2446 193770 4 1.717524.6953 0.4727 
Oe «43292 1.2429" 1.3719 4. 4582 mpd 2.6290 0.4715 
160 4.3212 1.2402 1.3672 4.4482 1.7138 2.6253 0.4736 
180 4.3160 1.2368 1.3662 4.4454 1.7130 2.6238 0.4762 
Puomeansi46 12333 123638 mas) “puseee e267 0.4842 
220 4.3127 1.2296 1.3595 iWan26 “FB GuenzNG2E9 0.4940 
240 4.3061 1.2257 1.3540 4.430" qw7eo4 256004 0.5007 
260% 4.2988 192220 113487 wets m7 2eN 2KeRTe 0.5061 
280 4.2905 1.2188 1.3444 4 Gag mam 273 «(nea 0.5085 
300 4.2825. «1.2165 1.3417 4.4077 1.7243 2.6188 0.5078 
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DISCUSSION OF RESULTS 


The elastic behavior of a cubic crystal is completely specified by 


the two shear moduli, C = C,, and C'= ) and the reciprocal com- 


1 
4b DRESS) 
pressibility or bulk modulus, B = 3 (C1420) 5)- 


sistance to shear deformation across the (100) plane in the [010] direct- 


C is a measure of the ree 


ion and C' is the resistance to shear deformation across the (110) plane 
in the [1i0] direction. B represents the resistance to deformation by a 
uniform hydrostatic pressure. 

For tantalum the computed values of C and C' are shown plotted against 
temperature as the solid lines of Figs. 6 and 7. Fig. 8 shows a corre- 
sponding solid line for B. These plots for C and C' show that these 
quantities vary with temperature in a manner similar to that observed for 
Other metals. They approach O°K with zero slope and depend linearly on 
temperature for a good portion of the temperature range. However B as a 
function of temperature (Fig. 8) shows a departure from previously ob- 
served behavior by displaying a positive slope up to about 100 K, then a 
slope reversal and negative slope up to room temperature. 

In obtaining the elastic constants from a perturbation calculation it 
would be expected that, if the anisotropy is not a strong function of tem- 
perature, the perturbation corrections would be small and of the same sign 
at all temperatures. Thus if the experimental points decrease monotonically 
with increasing temperature as in Fig. 1, the elastic constants which only 
differ by small amounts from the experimental data should exhibit the same 
temperature dependence. That this is not the case for B is shown by the 
solid line in Fig. 8 which is plotted from the elastic constants listed in 
Table II. It is possible that a fortuitous choice of propagation direction 
combined with the inevitable errors in angle and transit times may lead to 
the observed behavior. 

In order to observe the effects of errors of measurement or orientation 
on the approximation method used in obtaining the elastic constants of 
tantalum, the input data were slightly changed and the elastic constants 
recalculated using a CDC 1604 computer. Thus, the measured values of ev" 
were individually and collectively changed by + 1%, the polar angles being 


held constant, or a measured polar angle 9 or 9 was changed by + yas the 
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The calculated shear constant Cry, for tantalum 
perature (solid line). The dashed lines show C 


6, or @ are changed by the indicated amounts. 
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Fig. 7 The calculated shear constant $(C,4-C,5) for tantalum as a function of 


temperature (solid line). The dashed lines show $(c ) versus T 
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ie. 8 The calculated bulk modulus, $C, 426, 5), for tantalum as a function 
of temperature (solid line). The dashed lines show $(6, +26, 5) versus 


T when ev’, Q@, or @ are changed by the indicated amounts. 





values of ev" and the other angle being held constant. The results of 
these calculations are indicated as dashed lines in Figs. 6, 7, and 8 
and show that for substantial input data changes the curves for C, C', 
and B are merely shifted up or down their respective ordinates with no 
changes in shape. Because there is no change of shape with these per- 
turbations of the input data, it is concluded that the observed slope 
reversal effect in the bulk modulus curve for tantalum is a real effect. 

In the tungsten crystal C and C' were measured directly since pro- 
pagation was along [110] - These results are shown in Fig. 3. The 
variation of B with temperature has been calculated and is shown as the 
upper curve in Fig. 9. 

The most striking feature of the variations of the measured quanti- 


ties C, C' and $(C,,+C,,) +C/, with temperature is their very flat be= 


havior from liquid air down < liquid helium temperatures. This is the 
result of there being no observable difference in the elastic wave transit 
times for the three modes from liquid air to liquid helium temperatures. 
Because of this, the bulk modulus, B, remains constant with increasing 
temperature out to approximately 100°K. Beyond this point B displays a 
positive slope out to 160°K when it takes up a negative slope and a linear 
variation with temperature. 

Fig. 3 shows tungsten is nearly isotropic in the temperature region 
measured, being exactly so at the cross-over point between the two modes 
at about 195°K. Above 140°K the shear modes vary quite linearly with 
temperature. 

As was the case with tungsten, C and C' were measured directly for 
molybdenum. These results are displayed in Fig. 5, which shows that C 
varies monotonically with temperature but C' displays a tendency to level 
off in the 90°K and 180°K regions. 

The variation of B with temperature has been calculated and is shown 
as the lower curve in Fig. 9. It is seen that the slope of the bulk modulus 
curve goes from zero to negative as the temperature rises from o°K but there 
are two subsequent pairs of gradual slope reversals over the remainder of 
the temperature range measured, one reversal at about 80°K and one at about 


230°K. Reference to the 1% error flags in Figs. 5 and 9 readily establishes 
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that these slope variations of C' and B all lie within about 0.2% of a 
line that would show monotonically decreasing behavior with temperature. 
However as Fig. 4 indicates, two separate measurements of the combination 
$(C,,+C,,) +C) versus temperature (the longitudinal acoustic mode for 

[110] ) were made from liquid air to room temperatures, with consistent 
results. This consistency indicates that the depicted rather than the 
monotonic behavior of B versus T is warranted. 

It is seen from Fig. 1 - 5 that the temperature variation of the ob- 
served quantities is similar to that observed for other metals. That is 
the elastic constants approach 0°K with zero slope and, except for molybdenum, 
are remarkably linear over a wide temperature range. The percentage change 
of the constants relative to their 0°K values is, with the exception of 
tantalum, about the same as that observed for the noble metals. Near room 
temperature all the constants vary linearly with temperature and have 
Significantly smaller logarithmic derivatives than the noble metals. Un- 
fortunately, the change of elastic constants with pressure have not yet 
been measured for these materials so that it is not possible to evaluate 
the intrinsic temperature dependence of the elastic constants. 


> 


Other investigators have reported results for these metals obtained 


solely by dynamic methods. Also reported in the literature? > are re- 
sults obtained by dynamic measurements methods combined with Bridgman's 
compressibility dns, 2 

Table V compares room temperature values from the literature with the 
results of this paper. The differences between the tantalum results re- 
ported by Bolef and by this paper tend to be close to the sum of the re- 
spective estimates of error (~ 2%). This represents fair agreement. The 
tungsten results reported by Bolef and de Klerk are in quite good agree- 
ment with those reported in this paper (within 1%). 

We note that for tungsten the difference between Bridgman's value of 
the compressibility and the result of this work is very close to be cal- 
culated difference between the room temperature adiabatic and isothermal 


compressibilities. For tantalum and molybdenum, however, no such agree- 


ment is found, using either set of ultrasonically derived constants. 
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TABLE V. The reported room temperature values of the adiabatic elastic _ 
constants of tantalum, tungsten, and molybdenum, in units of 10 dyne-cm. 


Element @ (on B Investigator 
Tantalum 0.825 0.53 1.96 Bolef* 
0.818 0.518 1.919 This Paper 
wanesten i Sil Test 2.99 Were 
1.53 1.54 3.08 Bridgman’ 
1.604 1.598 3.083 Bolef & de ieee 
1.607 1.594 3.108 This Paper 
Molybdenum T.09 1.4 273 Druyvesteyn~ 
1.068 1540 2.683 Bolef & de ere 
226 1.342 2.619 This Paper 
a See Reference 8 
b See Reference 1l 
c See Reference 12 
d See Reference 9 


e See Reference 10 
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The difference in the values of the elastic shear constants reported 
for molybdenum by Bolef and de Klerk and by this paper are quite large. 

It seems clear, particularly in light of the agreement on tungsten, that 
these differences in moduli are real and that they depend on physical 
differences between the specimens rather than on differences in the mea- 
suring apparatus. That this must be the case is made more credible by 
considering that in measuring Chas which is measured directly for [110] 
propagation, the difference in the two reported values corresponds to a 
difference in transit time of 0.5 ps in the sample used, far outside the 
region of any reading error. 

A possible explanation of the difference in reported room temperature 
elastic constants (~12% differences for C and C' from the values of Bolef 
and de Klerk) for molybdenum is the influence of alloying of the impurities. 
Small amounts of a second phase have been found to change the single crystal 


a sive and 1#éfium°. These changes 


elastic constants of copper, 
in elastic constant have been interpreted in terms of changes in the long 
range and short range contributions to the elastic constants. The experi- 
ments indicate that the percent change in modulus per percent change in 
composition is seldom more than four and usually near unity. Since for 
molybdenum the observed percent change in modulus per percent change in 
composition is about 100, it is difficult to see how these differences can 
be interpreted on the basis of simple alloying. 

Anelastic effects due to dislocations also contribute to change the 
moduli. A dislocation contribution (diminution) to the moduli of about one 
percent has been observed for a single crystal of eonpeas: and of eae 
of purity in the 10 to 200 ppm range. In these experiments, the addition 
of small amounts of impurity tended to raise the observed modulus through 
dislocation pinning. Chambers and Schultz?” have observed an internal 
friction spectrum for molybdenum in the low temperature range that shows a 
pronounced peak in the region just below 100°K (the peak shifts with 
annealing) and a high level minimum at about 200°K. They explain certain 
of the observed internal friction effects as being due to the interactions 


of gaseous interstitials and dislocations. The assay of the present moly- 


bdenum sample (Table I) shows a significantly higher gas content than the 
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sample of Bolef and de Klerk. In addition, the present sample is much 

larger, and since dislocation densities increase rapidly with the dia- 

Meter of the speemen >. conditions are right in the present sample for 
an anelastic mechanism such as that proposed by Chambers and Schultz to 
be quite pronounced. 

The differences in the reported elastic shear constants for moly- 
bdenum are in opposite directions, and if these differences can be ex- 
plained by means of an anelastic mechanism, then these differences should 
yield some evidence as to how dislocations are distributed in the lattice. 
pooeer- has conducted an analysis of low temperature internal friction 
peaks wherein only a small fraction of the dislocations are considered to 
lie along crystallographic directions. The work of Thompson and Helmesee 
on copper indicates, however, that an appreciable fraction of the dislo- 
cations tend to lie along principal crystallographic directions. That 
the two orthogonal shear modes in the molybdenum sample are almost dia- 
metrically affected by this presumed dislocation-interstitial gas anelastic 
mechanism indicates that such a preferred alignment of dislocations may 


indeed be present. 
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BULK MODULUS 


As shown in Fig. 8 the bulk modulus of tantalum rises with increasing 
temperature above the 4.2°K value before beginning to decrease linearly. 
Observation of this anomalous behavior of B with temperature has led the 
authors to calculate the bulk moduli of the other transition elements 
whose elastic constants have been determined ultrasonically over the same 
temperature range. The following data were used: Alene for vanadium; 
Rayne and Ghandrasetnanas for iron; Alers, Neighbours, and aeons for 
nickel; Raynes’ for palladium; eacrolie! for niobium; and this paper for 
tantalum, tungsten, and molybdenum. Curves presenting these results are 
shown in Figs. 10 and 11 where the bulk modulus has been normalized to its 
0°K value, denoted by Bo: Fig. 10 shows B/B versus temperature for the 
body=centered cubic transition metals. Here the behavior of niobium, 
tantalum, tungsten, and molybdenum is significantly different from that 
displayed by vanadium and iron. Similarly in the face-centered cubic 
elements shown in Fig. 11 palladium has quite a different temperature de- 
pendence from that displayed by nickel. Raeaa has ascribed the departure 
from linearity at temperatures above 100°K of the elastic constants of 
palladium to the change with temperature of the contribution to the elastic 
constants from the unfilled d band. B/B | values from the elastic con- 
stants of apes ae silver and Pollde have been calculated and are included 
in Fig. 11 to indicate the temperature dependence of the bulk moduli of 
metals that have closed-d shells. 

Taking monotonically decreasing behavior with temperature as the norm, 
Figs. 10 and 11 show that the 3d and closed-d shell metals follow the norm. 
Niobium and palladium, 4d metals, change the sign of their slope several 
times over the temperature range measured, displaying a strong positive 
slope as room temperature is approached. Tantalum (5d) shows an early 
positive slope, then tends to follow the norm above 100°K. Tungsten (5d) 
has a flat response from O°K out to 100°K where a positive slope finally 
asserts itself. B/B for tungsten then parallels that of tantalum very 
closely on up to room temperature. Molybdenum (4d) generally follows the 
norm but with undulating departures. 

Figs. 10 and 11 indicate a significant difference between the behavior 


of the bulk moduli of the 3d metals, and the 4d-5d metals. Such a dif- 
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ference is in consonance with the observance of ferromagnetism in the 3d 
group but not the 4d-5d groups. 

One possible reason for this anomalous behavior of B with temperature 
that suggests itself is that these metals are involved in antiferromagnetic 
ordering in these temperature ranges, with associated anomalies in their 
thermal expansion coefficients. Evidence in support of this from other ex- 
periments is scarce and difficult to assess. Bolef, de Klerk, and Brandes: 
have reported anomalous thermal expansion behavior in the 180°K to 240°K 
region for vanadium. Fig. 10, however, shows that vanadium is a well-behaved 
3d metal in its B behavior. Hoare and Tenia o found a maximum in the 
Magnetic susceptibility of palladium between 80°K and 100°K yet Crangle and 
rien recently looked for a specific heat A -point in this region for 
Palladium, with no success. They also measured the specific heat of vana- 
dium, confirming the results of Clusius, Franzosini, and Picebetecnes that 
vanadium shows no anomaly between 175°K and 265°K, the region of Bolef, 
el al's anomalous thermal expansion behavior for the same metal. From the 
point of view that the body-centered cubic structure is the more natural 
vehicle for exhibiting antiferromagnetic behavior, and looking at the 
anomalous behavior of B/B versus T for the 4d-5d members of this group 
as possible evidence of this effect, then it is somewhat disconcerting that 
face-centered cubic palladium displays the same anomalous behavior. As the 
results of detailed specific heat, magnetic susceptibility, and thermal ex- 
Pansion experiments for pure samples of these metals become available, the 
role of antiferromagnetic ordering in the transition metals may become 
clearer. 

For body-centered and face-centered cubic structures de Penne has 


raised the possibility that C is equal to the bulk modulus of the 


127 °44 
electron gas of a metal. The temperature dependence of a free electron 
bulk modulus should show a steady decrease with increasing temperature as 
the lattice parameter increases. Reference to the Ci oP yy, columns in 
Tables II, III, and IV shows that tantalum, tungsten, and molybdenum like 
copper and silver exhibit just the opposite effect, an increase of Cio Suy 


with temperature. This reinforces de Launay's original admonition that the 


bulk modulus as derived from his considerations is a phenomenological quantity. 
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It does not reflect on the validity of the Sommerfeld free electron model. 
However correlation of Cio Sn, with the electronic bulk modulus of mono- 
valent metals is still considered a poeeibinic 

Excluded in the above calculations is the energy of the lowest state, 
Wo: Although its effect on the bulk modulus is usually taken to be small, 
it is possible that this may not be the case. Further consideration of 


this contribution await a cellular method calculation of Woe 
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OVERLAP FORCES 


35 : ; 
Fuchs has accounted for the elastic constants of copper as being 


composed of two main parts: the electrostatic stiffnesses due to that 
part of the energy which can be regarded as the energy of an array of 
positive point charges immersed in a uniform negative sea; and the over- 
lap energy due to the overlap of the charge distributions of ions plus 
the exchange energy. Following the lead of Fuchs and Znen Tsenhere: 
has concluded that the major contribution to the elastic constants of the 
body-centered cubic transition metals is due to the overlap energy (the 
Coulomb interaction energy due to the overlap of the charge distribution 
of ions plus the exchange energy due to the same overlap). 

The relations giving the overlap contributions to the elastic con- 


stants can be Piccen’ as: 


N [re w,'] . 


N [e“w,"! 1 3/4 (C, ,+2C 


3/8 (C,,-C, 9) 


4, 


rae 2 es 
N fr Woh oe Siiecn: 


where N is the density of atoms, W, is the nearest neighbor interaction 


1 


energy and W, is the next nearest neighbor interaction energy, and the 


primes Pain. differentiation. To a first approximation, the first 
derivative of Wo may be Bee creda a and the remaining derivatives are 
evaluated at the nearest of next-nearest neighbor distance as appropriate. 
In applying the above formulas, only the part of the elastic constants due 
to overlap forces is included, i.e., the overlap contribution to the modu- 
lus is taken to be the observed modulus minus the contributions from all 
other sources. 

Using the room temperature values of the elastic constants of tungsten 
as measured by Welene. Isenberg determined that the second derivative of 
the next-nearest neighbor interaction, Wei" > is positive, indicating a re- 
pulsive second-nearest neighbor interaction in tungsten. The values he 


obtained for the first and second derivatives of the nearest neighbor in- 


teraction Wy and Wo indicated a binding interaction for nearest neighbors. 
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Isenberg interpreted these results as a possible substantiation of Zener's 
>) of an antiferromagnetic arrangement for the body-centered cubic 
transition elements vanadium, chromium, niobium, molybdenum, tantalum, and 
tungsten. 

Taking the published values of the O0°K elastic constants of sauces 
mona. , and es and the results of this paper these same derivatives 
have been calculated. The results in units of ite dynes ae are shown 
in Table VI. For all these transition elements except iron the results are 
similar to Isenberg's. That is, the first derivative of the nearest neighbor 
interaction energy is negative and the second derivatives of both interaction 
energies are positive. 

The magnitudes of the second derivatives generally increase as one 
passes down and to the right in the periodic table. The magnitude of the 
first derivative of W increases to the right in the periodic table. The 
derivatives for iron are different in sign and magnitude as might be exe 
pected for a ferromagnetic material. 


1g 


; : , =< loa 
We note that if a simple exponential form W, = A.e i is chosen for 


both W, and W, then without neglecting W,' it can be shown 


Zi oz 
ce +Q, 4 ) va ibe C1 Ong 


HI 


9/8 (C, 5+C,,) 


[oye + ois =) | l 
where the subscript on the square bracket denotes the appropriate distance 
at which the quantity in the bracket is to be evaluated. The elastic con- 
stants used above are understood to include only that part due to overlap. 
For the transition metals considered here ore oe and C,,>0 so that the 
above relations predict Wy and Wo both to be repulsive energies. This re- 
sult is in agreement with ascribing the large cohesive energy of molybdenum 


: Boh be 
and tungsten to the large number of electrons in the low lying d band. : 
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TABLE VI. Values of the derivatives of the nearest and next-nearest 
neighbor interaction energies attributable to overlap forces, as 

fe) ' : 
calculated from the 0 K elastic constants. Here N = the density of 


atoms, W. = the nearest neighbor interaction energy, W, = the next- 


1 Z 

nearest neighbor interaction energy and the primes indicate differ- 
entiation with respect to r. The subscripts 1 and 2 outside the 
brackets indicate evaluation of the enclosed quantities at the 
nearest and next-nearest neighbor distances, respectively. The value 


given are in units of 1; dyne-cm 


Element N E "| 1 N My" 1 N [ev "| 9 
Vanadium -1.8 oe 15.4 
Niobium -3.2 TO.Z 202 
Tantalum -2.1 2072 i533 
Molybdenum -1.0 2602 PAS bsp 
Tungsten -0.8 34.5 3a 
Iron +0.8 20.2 pms 
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DEBYE TEMPERATURES 


Using Houston's method as outlined by Betts, Bhatia, and women the 


Debye temperatures of tantalum, tungsten, and molybdenum have been calcu- 
lated from their 0°K elastic constants. Table VII presents these results 
and compares them with the results of calorimetric measurements. Although 
the values entered for the Debye temperatures of tantalum in Table VII 
differ by rather large amounts, the lattice specific heat given by 

9 = 261.9°K is not in disagreement with the lower temperature part of the 
specific heat data of File et al and wore: et al. 

The data of Keesom and Deer eenee are too scattered to admit critical 
comparison and wisteone = lists only the value 8, = 245°K. Thus the 
situation for comparison of elastic and calorimetric Debye temperatures of 
tantalum issimilar to that for other etetet i.e., where the calorimetric 
data is not too scattered, the lattice specific heat calculated from the 
elastic constants fits the low temperature part of the calorimetric data. 

Much the same remarks apply to tungsten and molybdenum. For tungsten 
the data of Clusius and Reanzesinin is fitted at the low temperature end 
by oe = 384.3°K. Clusius and caus list 8. = 380°K and both Debye 
temperatures fit the data quite well. As in the case for tantalum Wolcott 
lists only a number and other data in the literature is too scattered for 
comparison. 

For molybdenum, the calculated value 8, Bahra oha'e gives a lattice 
specific heat in agreement with the lower temperature portion of the data 
of Clusius and Franzonini. The other published calorimetric data is not 


suitable for comparison. 
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TABLE VII. 


specific heat data. 





Element 8 5 Reference 
Tantalum 261. This Work 
252 a 
246. b 
245 c 
231" d 
Tungsten 384. This Work 
380 e 
405 c 
Molybdenum 474. This Work 
470 G 
454 e 
445 £ 
a See Reference 41 
b See Reference 43 
c See Reference 44 
d See Reference 42 
e See Reference 46 
f M. Horowitz and J.G. Daunt, Phys Rev 91, 1099 (1953) 
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The Debye temperatures of tantalum, tungsten, and molybdenum 
as calculated from the elastic constants and as derived from low temperature 
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